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Abstract 

An  analysis  of  mass  transfer  losses,  or  concentration  over-potentials  in  fuel  cells  is  provided.  An  elementary  theory,  based  on  an 
equivalent  film  thickness,  as  proposed  in  some  texts,  is  derived.  This  is  followed  by  a  more  rigorous  theoretical  treatment  of  mass  transfer 
theory,  for  which  the  mass  transfer  factor  is  obtained  as  a  function  of  the  driving  force.  The  solution  for  the  driving  force  is  derived,  for  the 
well-known  one-dimensional  convection-diffusion  problem.  It  is  shown  that  mass  transfer  in  planar  and  square  geometries  approximates 
this  idealised  situation.  A  linearised  theory,  appropriate  for  low  mass  flow  rates  is  also  presented.  The  methodology  is  illustrated  using  the 
simple  example  of  a  solid  oxide  fuel  cell  (SOFC).  It  is  shown  that  the  simplified  theory  is  only  applicable  for  very  dilute  binary  mixtures.  A 
step-by-step  procedure  for  computing  mass  transfer  in  fuel  cells  is  detailed,  together  with  a  discussion  of  the  scope  and  range  of  application 
of  the  results. 
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1.  Introduction 

When  current  is  drawn  from  a  fuel  cell,  the  external  volt¬ 
age,  V,  may  be  written  as 

V  =  E  -i"r  (1) 

where  E  is  the  ideal  or  Nernst  potential,  i"  the  current 
density* 1,  r  the  resistance,  and  rj  are  sometimes  referred  to 
as  ‘over-potentials’  or  ‘polarisations’. 

The  performance  of  fuel  cells  is  generally  accepted  to  be 
reduced  by  kinetics  (activation),  Ohmic  resistance  and  mass 
transfer  (concentration  effects).  At  low  current  density,  acti¬ 
vation  reduces  the  voltage,  at  intermediate  voltages,  Ohmic 
losses  within  the  electrolyte  are  dominant,  while  at  high  cur¬ 
rent  densities,  necessary  if  fuel  cells  are  to  generate  large 
sources  of  power,  mass  transfer  effects  predominate.  Fig.  1 
shows  a  current-voltage  distribution  curve  typical  of  many 
hydrogen  fuel  cells.  In  this  paper,  we  are  concerned  with 
the  mass  transfer  effects,  which  predominate  at  high  current 
densities.  While  numerous  texts  on  fuel  cells  are  available 
today,  many  of  these  are  somewhat  vague  on  the  mathemat¬ 
ical  details  needed  to  calculate  mass  transfer  effects.  The 
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1  In  this  manuscript,  the  convention  of  Jacob  [1]  whereby  a  ‘dot’ 
represents  a  time  derivative,  and  a  ‘dash’  a  space  derivative,  is  adopted. 


purpose  of  this  paper  is  to  provide  the  reader  with  a  reason¬ 
ably  simple  yet  comprehensive  method  by  which  he  or  she 
may  estimate  these  factors. 

Historically,  chemists  and  chemical  engineers,  have  con¬ 
sidered  mass  transfer  problems  based  on  molar  concentra¬ 
tions,  Xi,  see  for  example  the  book  by  Bird  et  al.  [2].  Using 
this  latter  convention,  diffusion  losses  may  be  introduced 
into  Eq.  (1)  in  the  generic  form 


RT 

nvF 


In 


where  +b  and  +w  are  the  molar  concentrations  in  the  bulk  of 
the  fluid  and  near  the  wall,  n  an  integer  which  is  a  function 
of  the  stoichiometry  of  the  reaction(s),  v  the  valence  of  the 
ionic  conductor,  and  F  is  the  Faraday’s  constant. 

The  theoretical  foundation  for  the  present  manuscript 
follows  the  work  of  Spalding  [3].  This  is  based  on  mass 
fraction,  ra;  and  mass  flux  relative  to  mass-average  veloc¬ 
ity,  j/' .  A  justification  for  the  use  of  these  widely  adopted 
definitions  is  that  in  the  conservation  laws  of  continuum  me¬ 
chanics,  it  is  mass  (rather  than  mole  numbers)  momentum 
and  energy  which  are  conserved,  and  indeed  mass-based 
formulations  and  methodologies  and  codes  abound  in  the 
engineering  community.  Mass-based  and  mole-based  ap¬ 
proaches  are  easily  inter-converted  as  follows: 

Mi 

mi  =  Xi —  (3) 

M 
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where  the  mixture  molecular  weight  is  just 

1 


M  =  J2  X,  M,  = 


YJmi/Mi) 


(4) 


In  the  present  formulation,  the  diffusion  losses  may,  thus 
be  expressed  as 


RT  /mb 

rj  =  - In  I  - 

nvF  \mw 

for  reactants,  and 


RT 

ii  =  - In 

nvF 


mh 


(5) 

(6) 


for  products,  where  m b  and  mw  are  bulk  and  wall  values 
for  mass  fraction.  The  purpose  of  this  paper  is  to  derive  ex¬ 
pressions  which  allow  for  Eqs.  (5)  and  (6)  to  be  expressed 
analytically.  The  terms  ‘over-potential’  or  ‘polarisation’  are 
widely  used  and  not  entirely  unambiguous.  Moreover,  the 
term  RT/nvF ,  serves  very  little  purpose  for  a  fundamental 
analysis.  Therefore,  we  shall  define  a  mass  transfer  factor  as 
ln(mb/mw),  for  reactants,  and  ln(mw/mb),  for  products.  De¬ 
termination  of  the  mass  transfer  factors  for  the  components 
of  the  fuel  and  air  mixtures,  is  the  goal  of  this  work. 


m"  =  ± 


Faraday’s  law  may  be  written  as 
Mi" 


vF 


(10) 


i”  =  ±—g(mv/  -  mb) 


The  current  density  is,  therefore,  given  by 
Fv 


(11) 


For  wall  suction,  the  maximum  possible  current  density 
is  obtained  when  mw  — >  0 

.//  Fv 

i  max  —  (12) 


Elementary  mass  transfer  analyses  employ  an  ‘equivalent- 
film-theory’  to  obtain 


(13) 


where  8  is  the  thickness  (length- scale)  of  an  imaginary  film. 
Since  by  Eq.  (13),  the  conductance,  g,  is  constant 


mb 


r' 

1  max 


(14) 


and  Eq.  (7)  is  obtained,  identically.  We  shall  refer  to  Eq.  (7) 
as  the  ‘simplified  approach’  in  the  remainder  of  the  text. 

This  analysis  is  not  recommended  by  the  present  writer 
for  a  number  reasons:  (i)  for  high  mass  transfer/current 
densities,  the  conductance,  g,  may  be  a  function  of  /",  and 
not  constant;  (ii)  Fick’s  law  has  been  written  in  terms  of 
absolute  rather  than  relative  velocity;  (iii)  the  notion  of  a 
conductance  is  generally  introduced  in  a  form  such  that 
the  diffusion  flux,  j"  =  F(dm/dy),  may  be  replaced  by 
g(mw  —  mb),  an  expression  for  the  mass  flux,  m"  =  pu , 
requires  the  introduction  a  mass  transfer  driving  force,  or 
equivalent.  This  is  discussed  further  below.  On  a  more  prac¬ 
tical  note,  the  equivalent  film  model  does  not,  in  any  way, 
correspond  to  the  physical  reality  of  flow  and  mass  transfer 
within  the  passages  of  a  fuel  cell,  which  are  typically  ducts 
of  rectangular,  or  other,  cross-section. 


1.1.  Simplified  analysis  based  on  equivalent  film  theory 


Some  authors  [4,5]  have  argued  that  diffusion  effects  in 
fuel  cells  may  be  modelled  using  terms  of  the  form: 


RT  (  i"  \ 

- In  1  -  - - 

nvF  \  /"max  / 


where  /"max  is  a  ‘limiting  current  density’. 

A  typical  argument  proceeds  as  follows  [6]:  the  rate  of 
mass  transfer  may  be  represented  by  Fick’s  law 

.  „  dm 

m"  =  pu  =  F —  (8) 

dy 

where  r  =  pD  is  an  exchange  coefficient.  It  is  generally 
useful  to  introduce  the  concept  of  a  mass  transfer  conduc¬ 
tance,  g,  by  the  following  rate  equation: 

m"  =  g(mw  -  mh)  (9) 


2.  Elements  of  mass  transfer  theory 

2.1.  Definition  of  the  driving  force  and  blowing 
parameter 

The  present  section  follows  that  in  [3,7],  the  latter  of 
which  contain  reviews  of  Bird  et  al.  [2],  Treybal  [8],  Sher¬ 
wood  et  al.  [9],  and  other  works.  The  methodology  has  been 
advanced  in  a  number  of  more  recent  texts  [10,11].  In  this 
paper,  the  theory  is  further  refined  for  application  to  fuel 
cells. 

In  fuel  cells,  heterogeneous  chemical  reactions  occur  on 
electrode  surfaces.  These  lead  to  sources  and  sinks  in  the 
continuity  and  species  (mass  fraction)  equations.  The  mass 
flux,  m",  is  given  by 

m"  =  gB 


(15) 


S.B.  Beale/ Journal  of  Power  Sources  128  (2004)  185-192 


187 


Bulk  state 


m, 


rom 

mm  I  — 

W 

oy 


Wall 


t 


T . t 


w 


m. 


Control 

volume 


T- state 


t 


mx 


rii'm 


Fig.  2.  Notion  of  ‘transferred  substance’  state. 


where  B  is  a  mass  transfer  driving  force 

mb  -  mw 

li  =  - 

mw  —  mt 


(16) 


Both  m"  and  B  are  positive  for  injection  and  negative  for 
suction,  while  the  conductance  g  is  always  positive.  The 
subscripts  refer  to  values  in  the  bulk,  at  the  wall  and  at 
the  so-called  4 transferred- substance  state’  (T-state).  Fig.  2 
illustrates  the  T-state  concept,  schematically.  For  multi- 
component  mixtures,  T-state  values  may  be  computed  as 


(17) 


where 

m"  =  m"  j  (18) 

j 


In  the  absence  of  chemical  reactions,  for  example  mass 
transfer  in  reverse  osmosis  membranes,  the  T-state  has  phys¬ 
ical  significance,  namely,  a  state  far  enough  from  the  wall 
such  that  diffusive  effects  are  negligibly  small,  i.e.  mt  is  a 
reference  mass  fraction,  and  0  <  mt  <  1.  For  situations  such 
as  the  present,  where  chemical  reactions  are  present,  the  no¬ 
tion  of  the  T-state  and  mt  is  somewhat  abstract,  though  no 
less  useful,  with  —  oo  <  m{  <  ±oo. 

Eq.  (15)  reduces  to  Eq.  (9)  for  two  limiting  cases  [7]: 

(i)  5  — >  0,  corresponding  to  low  mass  transfer  rates  asso¬ 
ciated  with  very  dilute  mixtures. 

(ii)  mt  — >►  ±oo,  this  situation  occurs,  for  example  in  chemi¬ 
cal  catalytic  reactors,  where  there  is  no  net  mass  transfer 
at  the  wall,  although,  there  is  still  internal  diffusion. 


Under  all  other  circumstances,  it  is  Eq.  (15)  which  should 
be  used,  and  Eq.  (9)  avoided.  It  is  perhaps  unfortunate  that 
the  term  'low  mass  transfer  theory’  has  evolved  in  the  lit¬ 
erature  since,  as  will  be  shown  below,  there  are  situations 
where  mass  transfer  rates  are  small,  which  do  not  correspond 
to  cases  (i)— (ii).  Conversely  for  case  (ii),  there  can  still  be 
large  diffusion  fluxes  and  gradients. 

For  convenience  let  a  blowing  factor,  b ,  be  defined  by 


(19) 


where  g*  is  the  value  of  the  conductance  in  the  limit  m"  ->  0 
(NB:  g*  is  often  obtained  from  heat  transfer  analysis).  Com¬ 
bining  Eq.  (19)  with  Eq.  (15),  we  obtain 

B  =  —b  (20) 

8 


which  is  the  non-dimensional  ‘Ohm’s  law’  of  mass  trans¬ 
fer,  namely  that  the  driving  force,  B ,  is  proportional  to 
the  blowing  parameter  (normalised  convection  flux),  with 
the  constant,  g*/g,  representing  the  normalised  resistance. 
Characterisation  of  any  two  of  the  three  non-dimensional 
numbers  in  Eq.  (20)  is  a  sufficient  solution  of  the  mass 
transfer  problem.  In  practice,  Eq.  (20)  is  non-linear  since  g 
is  a  function  of  b. 

It  is  simple  to  rearrange  the  definition  of  the  driving  force 
Eq.  (16),  to  obtain  an  expression  of  the  required  form  for 
Eqs.  (5)  and  (6).  For  reactants: 


RT 

- In 

nvF 


/ 1  +  r/A 

vi+U 


(21) 


For  products, 


RT 

- In 

nvF 


1  +  B 
1  +  rB 


(22) 


where 

mt 

r  =  — 
mb 


(23) 


In  conventional  mass  transfer  problems  it  is  g/g*  as  a 
function  of  B  (or  b)  which  is  specified.  What  is  required  for 
the  present  analysis  is  data  of  B  as  a  function  of  b.  This 
can  be  obtained  from  experiment,  by  detailed  numerical  cal¬ 
culations,  or  from  theoretical  analysis,  for  certain  idealised 
cases.  We  shall  consider  the  latter,  and  show  that  it  is  a 
reasonable  solution,  under  many  circumstances,  for  actual 
fuel-cell-duct  geometries. 


2.2.  One- dimensional  convection  diffusion 


One-dimensional  (ID)  convection-diffusion,  sometimes 
referred  to  as  a  Couette-flow  or  falling-film  problem,  may 
be  posed  in  the  form 

.  //  dm 

m  mt  =  pum  —  F —  (24) 

d y 

where  m  =  m(y).  This  equation  may  readily  be  solved 
using  the  method  of  separation  of  variables.  The  result  in 
non-dimensional  form  is 

B  =  exp  (b)  —  1  (25) 

Fig.  3  shows  B(b)  computed  using  Eq.  (25).  Also  shown 
are  numerical  solutions  to  the  convection  diffusion  problem 
[12]  for  fully-developed  fluid  flow  and  mass  transfer  in  plane 
and  square  ducts  corresponding  to  fuel  cell  channels  away 
from  the  entrance  region.  It  can  be  seen  that  the  results  of 
the  calculations  agree  well  with  the  ID  theoretical  curve. 
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b  =  m"/g* 


Fig.  3.  Driving  force  vs.  normalised  mass  flux. 


The  compression  of  these  data  onto  a  single  curve  is  due  to 
the  choice  of  non-dimensional  parameter(s),  and  would  not 
occur  if  the  abscissa  were  non-dimensionalised  in  the  form 
of  a  Sherwood  number 

gDh 

Sh=^jf-  (26) 

where  D h  is  a  hydraulic  diameter  [10].  The  reader  will 
appreciate  that  although,  the  results  of  the  numerical  data 
are  in  reasonably  close  agreement  with  Eq.  (25),  val¬ 
ues  of  Sh*  are  quite  different  for  the  three  cases;  with 
Sh*  equal  to  (a)  8.23,  (b)5.38,  (c)2.84  [13].  Fig.  4  shows 
values  of  Sh*  for  rectangular  ducts,  2 a  x  2 b,  typical  of  many 
fuel  cell  designs  as  a  function  of  aspect  ratio,  a  =  a /h.  This 
was  obtained  from  Table  44  in  the  book  by  Shah  and  Lon¬ 
don  [13].  Eq.  (25)  combined  with  Eqs.  (21)  and  (22)  is  the 
recommended  method  for  computing  mass  transfer  effects 
in  fuel  cells.  It  is  referred  to  below  as  the  ‘present  method’. 

2.3.  Linearised  method  for  low  mass  flux 

Fig.  3  also  shows  the  low  mass  flow-rate  solution,  i.e. 
\im\m\"  0,  for  which  it  may  readily  be  shown  that 

4  =  1  (27> 


Fig.  4.  Sh *  for  a  rectangular  duct.  Adapted  from  Shah  and  London  [13]. 


B  =  b  (28) 

This  is  referred  to  below  as  a  ‘linearised  solution’,  and 
may  only  be  applied  to  situations  for  which  K<  1,  a  situa¬ 
tion  which  occurs  in  many  fuel  cells  operating  at  present-day 
current  densities. 


3.  Example:  a  solid  oxide  fuel  cell  with  H2  as  fuel 

3.1.  Problem 


The  proposed  methodology  will  be  presented  by  means 
of  a  quantitative  example.  The  simple  case  of  a  solid  oxide 
fuel  cell  (SOFC),  with  H2  as  fuel  is  considered.  Oxidation 
takes  place  at  the  anode  surface: 

H2  +  O2-  -*  H20  +  2e“  (29) 

and  reduction  takes  place  at  the  cathode  surface: 

02  +  4e~  -*  202~  (30) 


The  Nernst  equation  for  a  solid  oxide  fuel  cell  may  be 
written  in  terms  of  molar  fractions  as 


E  =  E°  + 


RT 

IF 


In 


*h2*q25 
-xh2o 


RT 

H - In  Ea 

4  F 


(31) 


Thus,  v  =  2  (for  O2),  and  n  —  1,  2,  1  for  H2,  O2,  H2O, 
respectively.  The  Nernst  equation  may  easily  be  converted 
to  a  mass  fraction  form: 


E  =  E°  + 


RT 
2 F 


In 


RT 

H - In  Pa 

4  F 


'«h,'«o25 
mH2o 


where 


1 

-  In (Ma)  -  0.4643 
2 


(32) 


/  Mh2Mq5\  1 

In  - &  In (2)  +  -  ln(32)  -  ln(18)  =  -0.4643 

\  MH2o  )  2 

(33) 

Assume  the  oxidant  to  be  a  mixture  of  O2  and  atmospheric 
nitrogen.  Since  only  O2  is  transferred  at  the  cathode,  rat  =  1 , 
Fig.  5.  Fig.  6  shows  the  corresponding  situation  at  the  anode, 
if  the  fuel  is  pure  hydrogen,  the  reaction  may  be  treated  as 
a  ‘simple  chemical  reaction’ 

H2  +  O2-  -*  H20  +  2e-  n4 

1  kg  +  rkg  ->  (1  +  r)kg 

If  the  weight  of  the  electrons  is  neglected,  r  =  8.  Thus,  from 
Eq.  (17)  mt  =  —1/8  for  H2,  and  for  H2O,  mt  =  9/8.  On 
the  anode  side,  for  H2,  Eq.  (21)  applies  whereas  for  H2O  it 
is  Eq.  (22)  which  must  be  used. 
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Fig.  5.  Mass  transfer  at  cathode. 
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Fig.  6.  Mass  transfer  at  anode. 


3.2.  Solution 

Figs.  7  and  8  are  plots  of  the  mass  transfer  factor  versus 
blowing  parameter,  b ,  for  O2  on  the  cathode  side.  Fig.  7 
is  a  comparison  of  the  present  and  simplified  methods,  for 
suction  b  is  negative  and  for  O2  m b  >  mw  and  raw  >  rat.  The 
simplified  method,  Eq.  (7)  has  been  plotted  in  a  functionally 


f 


In 


V 


Fig.  7.  Present  method  compared  with  simplified  method  for  oxygen 
(cathode). 


Fig.  9.  Present  method  compared  with  simplified  method  for  hydrogen 
(anode). 


similar  form  by  substituting  for  the  independent  variable 
according  to 


b 


m 


a 


:// 


b 


max 


/fi//max 


V 

1  max 


(35) 


where  /?max  =  —  gm^  ( g  occurs  in  both  numerator  and  de¬ 
nominator  of  Eq.  (35),  and  is  assumed  constant  in  Eq.  (13) 
for  the  equivalent  film  thickness  theory;  therefore  it  does 
not  enter  the  calculations).  Fig.  8  shows  a  comparison  of 
the  present  analysis  with  that  obtained  using  the  linearised 
method. 

Figs.  9-11  show  anode  side  results.  Fig.  9  shows  a  com¬ 
parison  of  the  simplified  and  present  methods  for  H2  on 
the  anode  side,  for  the  fuel-side  there  is  blowing  and  b  is 
positive,  with  mt  >  mw  >  m b  for  H2.  Fig.  10  shows  a  com¬ 
parison  for  H2  between  the  present  and  linearised  analyses. 
Fig.  11  shows  a  similar  comparison  for  H2O. 


3.3.  Discussion 


The  results  of  Fig.  7  suggest  that  the  simplified  approach 
agrees  with  the  present  procedure  only  when  is  very 
small.  Inspection  of  Fig.  7  reveals  that  unless  the  O2  mass 
fraction  is  very  small,  the  simplified  methodology  is  not  in 
agreement  with  the  present  method,  even  at  very  low  values 


Fig.  8.  Present  method  compared  with  linearised  method  for  oxygen 
(cathode). 


Fig.  10.  Present  method  compared  with  linearised  method  for  hydrogen 
(anode). 
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Fig.  11.  Present  method  compared  with  linearised  method  for  water 
(anode). 

of  the  mass  flux,  fn ",  corresponding  to  low  current  densities. 
Thus,  this  approach  would  appear  only  to  be  beneficial  if  the 
air  were  highly  oxygen-lean,  a  situation  which  is  unlikely 
to  be  found  in  modern  fuel  cells,  which  typically  employ 
atmospheric  or  enriched  air.  Of  course,  it  is  true  that  the 
conductance  may  be  redefined  as  g  =  g/(mt  —  raw),  so  that 
rh"  =  g/(mw  —  mb).  That  however,  does  not  fundamentally 
change  the  argument. 

Larminie  and  Dicks  [4]  noted  problems  with  the  simplified 
analysis,  and  recommended  the  use  of  an  empirical  relation 
of  the  form 

rj  =  A\n(l-  -jf-)  (36) 

where  values  of  the  constant  A  differ  from  RT/nvF  by  al¬ 
most  an  order  of  magnitude.  There  should  be  no  need  to 
employ  empirical  relations,  the  existing  theory  is  sufficient 
to  adequately  describe  the  situation.  For  H2  on  the  anode 
side,  Fig.  9,  the  effect  of  the  T-state  value  not  being  unity  is 
such  that  the  simplified  method  is  in  complete  disagreement 
with  the  present  approach,  even  when  m b  is  very  small. 

The  results  of  Figs.  8  and  10  do,  however,  suggest  that 
there  is  good  agreement  between  the  present  method  and  the 
linearised  method  for  all  values  of  mb,  provided  the  mass 
flux  fn "  (i.e.  b)  is  sufficiently  small,  it  can  be  seen  that  the 
two  sets  of  O2  curves  are  in  agreement  as  b  — >►  0,  for  large 
b  there  are  significant  departures.  The  trend  is  repeated  in 
Fig.  1 1,  for  H2O.  Water  unlike  H2  and  O2  is  being  produced 
rather  than  consumed  by  the  reaction,  and  hence  there  is  no 
limiting  current/mass  flux.  Rather  the  build-up  of  water  has 
decreasing  impact  of  the  mass  transfer-limiting  effects  on 
the  analysis.  Of  course,  this  does  not  mean  the  production 
of  H2O  is  limitless,  since  the  current  density  is  limited  by 
the  mass  transfer  of  H2  and  O2.  When  computing  the  cell 
potential  in  Eq.  (1)  the  sum  of  the  three  contributions,  H2, 
O2  and  H2O  must  be  considered. 

For  many  fuel  cells,  the  linearised  approach  will  suffice, 
but  this  will  not  be  the  case  as  the  limiting  current  (mass  flux) 
is  approached  and  rh"  increases.  Thus,  the  linearised  analysis 
will  not  be  applicable  at  high  current  densities,  and  the  user 


should  always  perform  sample  calculations  to  obtain  trial 
values  of  b  prior  to  using  the  linearised  method.  The  present 
analysis  generates  meaningful  results  at  all  mass  transfer 
rates.  It  is,  therefore,  recommended.  While  it  is  possible  to 
construct  more  complex  problems,  this  simple  example  has 
served  to  illustrate  the  important  principles  of  the  analysis. 


4.  Method  for  calculating  rate  of  mass  transfer  in  fuel 
cells 

The  analysis  can  readily  be  used  to  estimate  mass  transfer 
in  fuel  cells,  at  all  mass  flow  rates  (current  densities),  low 
or  high.  The  following  approach  is  suggested. 

(1)  Estimate  g*  for  both  the  anode  and  cathode  geometries 
under  consideration  (this  is  required  to  compute  b  for 
both  air  and  fuel  sides).  For  many  duct  geometries,  un¬ 
der  conditions  of  fully-developed  flow,  these  are  widely 
published  in  the  form  of  zero  mass  transfer  Sherwood 
(Nusselt)  numbers,  see  for  example  Fig.  4.  Hence,  com¬ 
pute  g*  and  b  =  Mi" /vFg*. 

(2)  Compute  the  driving  force  B(b).  Best  practice  would 
be  for  this  to  be  obtained  from  measured  experimental 
data,  or  detailed  numerical  calculations  for  the  particular 
geometry/flow  conditions  under  consideration.  If  this  is 
not  available,  Eq.  (25),  should  be  used  as  an  estimate. 

(3)  Compute  the  T-state  value  from  the  stoichiometry,  and 
hence  compute  the  mass  transfer  factors  (i.e.  polarisa¬ 
tions,  77)  given  the  bulk  value,  using  Eq.  (21)  or  (22). 

(4)  Repeat  the  above  for  all  participating  species  for  the 
given  working  fluid  (fuel  or  air)  and  sum  the  terms  as 
per  Eq.  (1). 


5.  General  discussion 

The  mass-based  conductance,  g,  is  by  and  large,  equiv¬ 
alent  to  the  molar  conductance  the  symbol,  k*  =  g/M ,  in 
Bird  et  al.  [2]  (the  radical  is  to  remind  the  reader  that  k  is 
not  constant).  Reference  [7]  notes  the  mole  numbers  are  not 
preserved  in  many  chemical  reactions,  and  suggests  this  as  a 
reason  for  the  use  of  a  mass-based  formulation.  The  present 
author’s  position  is  that  provided  the  mixture-molecular 
weight  at  the  wall  and  in  the  bulk  are  nominally  equal, 
Mw  =  Mb,  molar-based  and  mass-based  formulations 
should  be  entirely  equivalent.  If,  however,  Mw  ^  Mb,  the 
issue  of  replacing  local  values  for  the  diffusion  fluxes  by 
global  (i.e.  integrated)  values,  g(mw  —  m\f)  and  k*(xw  —  Vb), 
may  not  render  identical  results. 

In  this  paper,  the  example  of  a  SOFC  with  H2  as  fuel  was 
given.  The  present  method  may  readily  be  extended  to  other 
fuel  cells,  involving  different  fuels,  more  complex  heteroge¬ 
neous  and/or  homogeneous  reactions,  thermal  radiation,  etc. 
Kays  and  Crawford  [10]  suggest  that  the  present  analysis  is 
appropriate  for  non-dilute,  non-binary  mixtures  provided  all 
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diffusion  coefficients  are  equal 

A2  =  A3  =  ...  Ty  =  T  =  pD  (37) 

This  situation  would  appear  to  be  reasonably  true  for  the 
high-temperature  gases  in  SOFCs.  It  may  also  be  approxi¬ 
mately  true  for  lower  temperature  PEMFC’s.  However,  the 
analysis  is  not  valid  for  liquid-based  electrolytic  solutions 
and  other  possible  scenarios.  Wilke  [15]  provides  an  approx¬ 
imate  method  for  computing  A7  f°r  non-binary  mixtures, 
for  a  given  driving  force,  which  is  exact  for  dilute  mixtures, 
and  has  enjoyed  widespread  use. 

It  might  be  argued,  that  there  are  situations  for  which 
the  present  analysis  does  not  hold,  and  for  which  there  is 
no  other  recourse  than  to  solve  the  Stefan-Maxwell  equa¬ 
tions  for  the  multi-component  system.  Knuth  [14]  considers 
circumstances  under  which  it  is  prudent  to  substitute  ordi¬ 
nary  diffusion  for  the  Stefan-Maxwell  system  of  equations 
for  non-binary  non-dilute  systems  where  A/  7^  A k-  Coffee 
and  Heimerl  [16]  considered  five  approximate  solutions  to 
the  Stefan-Maxwell  equations,  for  laminar  flames,  includ¬ 
ing  [15]  and  noted  surprisingly  minor  differences  in  com¬ 
puted  values.  Amali  et  al.  [17]  considered  Stefan-Maxwell 
and  Fickean  diffusion  in  porous  media  and  observed  differ¬ 
ences  of  4-10%.  Thus,  while  the  use  of  Fick’s  law,  and  the 
important  additional  step  of  replacing  the  gradient  diffusion 
term  with  a  rate  term,  may  not  necessarily  be  ‘correct’  in  the 
strict  sense  for  non-dilute  multi-component  mixtures,  it  does 
offer  a  reasonable  balance  between  mathematical  rigour,  and 
computing  effort. 

There  are  also  a  number  of  other  important  factors  to 
consider  when  analysing  mass  transfer  in  fuel  cells.  For  ex¬ 
ample  the  present  approach  must  be  modified  if  large  prop¬ 
erty  variations  occur;  for  example  the  Schmidt  number  may 
vary  as  a  result  of  temperature  gradients  in  the  passages 
of  the  fuel  cell.  Surface  diffusion  and  Knudsen  diffusion 
may  be  present,  and  ‘ordinary’  diffusion  is  not  sufficient 
to  characterise  these  processes.  While  the  importance  of 
Stefan-Maxwell,  surface,  and  Knudsen  diffusion  is  not  to  be 
understated,  this  should  not  be  used  as  an  excuse  to  justify 
the  use  of  over-simplified  methods,  of  questionable  value 
in  place  of  the  present  approach.  There  is  clearly  range  of 
applications,  where  the  present  method  will  give  a  good  es¬ 
timate  of  behaviour.  The  goal  of  this  work  was  to  provide  a 
reasonably  simple,  yet  versatile  methodology  by  which  the 
fuel  cell  designer  could  estimate  the  impact  of  mass  transfer 
on  their  design. 

In  this  paper,  the  focus  has  been  on  the  analysis  of  mass 
transfer  in  gas-flow  channels  of  fuel  cells.  Many  fuel  cells 
also  contain  gas  diffusion  and  electrode  layers  in  the  form 
of  consolidated  porous  media.  Many  methods  consider  dif¬ 
fusion  by  defining  an  effective  diffusion  coefficient,  with 

Aft  =  r-  (38) 

r 

where  s  is  the  void  fraction  and  r  is  a  tortuosity  or 
length-scale  (normalised).  While,  such  an  approach  may 


well  approximate  within-phase  diffusion  in  gas  diffusion 
layers;  in  electrode  layers  it  is  also  necessary  to  compute 
bulk-to-wall  mass  transfer,  by  calculating  a  driving  force, 
B ,  as  described  above,  since  mass  transfer  can  only  occur 
in  the  presence  of  a  concentration  gradient  between  the 
bulk  of  the  fluid  and  the  wall.  Thus,  it  is  again  necessary  to 
obtain  S7z*  (or  Nu *),  for  the  porous  media  and  also,  if  pos¬ 
sible,  the  B  versus  b  characteristic  curve,  prior  to  applying 
the  methodology  above. 


6.  Conclusions 

An  analysis  based  on  the  well-established  concepts  of 
mass  transfer  driving  force,  B ,  and  transferred  substance 
state  was  adapted  to  provide  a  method  for  estimating  the 
mass  transfer  factors  (concentration  polarisations)  in  fuel 
cells.  The  method  is  based  on  knowledge  of  B  as  a  function 
of  the  non-dimensional  mass  flux  (i.e.  current  density),  b , 
together  with  data  for  the  zero  mass  transfer  conductance,  g*, 
obtained  from  an  appropriate  Sherwood  number  correlation, 
Sh*,  for  the  geometry  under  consideration. 

It  was  demonstrated  that  the  use  of  a  simple  equivalent 
falling-film  method,  to  estimate  mass  transfer  is  unreliable, 
even  for  situations  where  the  rate  of  mass  transfer  is  very 
small.  A  linearised  analysis  does  provide  meaningful  results 
for  small  current  densities,  corresponding  to  \b\  <<C  1,  al¬ 
though,  it  does  not  appear  to  offer  any  major  advantages 
over  the  present  analysis,  which  is  valid  for  both  low  and 
high  mass  flow  rates.  A  methodology  whereby  calculations 
on  mass  transfer  factors  in  practical  fuel  cell  geometries  may 
be  performed.  This  has  the  advantage  of  being  in  a  form 
compatible  with  methods  used  elsewhere  in  the  analysis  of 
transport  phenomena  in  engineering  equipment. 
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